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Abstract — this paper presented comparison of the time specification performance between two type of controller for a 
Double Inverted Pendulum system. Double Inverted Pendulum is a non-linear .unstable and fast reaction system. DIP is 
stable when its two pendulums allocated in vertically position and have no oscillation and movement and also inserting 
force should be zero. The objective is to determine the control strategy that to delivers better performance with respect to 
n angle's and cart position. In this paper simple multi PD controller designed on the theory of pole placement 
and its performance is compared with Linear Quadratic Regulator controller using MATLAB and Simulink. 
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I. INTRODUCTION 

Te inverted pendulum offers a very good example for control engineers to verify a modern control theory. This can 
be explained by the facts that inverted pendulum is marginally stable, in control sense, has distinctive time variant 
mathematical model. The double inverted pendulum is a highly nonlinear and open-loop unstable system. The inverted 
pendulum system usually used to test the effect of the control policy, and it is also an ideal experimental instrument in the 
study of control theory [1, 2]. To stabilize a double inverted pendulum is not only a challenging problem but also a useful 
way to show the power of the control method (PID controller, neural network, PLC, genetics algorithm, etc.). 

In this paper common control approaches such as the linear quadratic controller (LQR) and PD controller based on 
a pole placement technique to overcome the problem of this system require a good knowledge of the system and accurate 
tuning obtain good performance [3-5]. This paper presents investigations of performance comparison between modern 
control and PD control for a double inverted pendulum system. Performance of both controller strategies with respect to 
pendulums angle and cart position is examined. 

II. MODELING OF DOUBLE INVERTED PENDULUM 

To control this system, its dynamic behavior must be analyzed first. The dynamic behavior is the changing rate of 
the status and position of the double inverted pendulum proportionate to the force applied. This relationship can be explained 
using a series of differential equations called the motion equations ruling over the pendulum response to the applied force. 
The double inverted pendulum is shown in Fig 1. The meanings and values of the parameters for inverted pendulum are 
given in Table 1. [4] 




Figl: schematic diagram of Double Inverted Pendulum 

To derive its equations of motion, one of the possible ways is to use Lagrange equations [6] 

d dL dl 
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table i : Parameters of double inverted pendulum 



Mass of thecart,( first pole, 
5.8kg(1.5kg,.5kg,.75kg) 



L1G1), LjOz) Length of pendulum first(21 1 ) and length of 

second pendulum ( ?1 : ) .. lin. L5m 

_g Center of gravity 9,8m/s 

F | Force applied to cart 



Where L = T - ^ 
coordinates q and n 
of the system are gi 

T= i(m + " 
2m 3 l 1 )xti 1 cos6 



' is a Lagrangian, Q is a vector of generalized forces (or moments) acting in the direction of generalized 
)t accounted for in formulation of kinetic energy T and potential energy V. Kinetic and potential energies 
;iven by the sum of energies of cart and pendulums. 

i 1 +m 2 +m 3 )x 2 + (^mjj 2 +2m 2 l 1 2 + 2m 3 l 1 2 )8 1 2 +^m 2 l 2 2 8 2 2 + (m 1 ! 1 +2m 2 ! 1 + 
\ + m 2 l 2 xti 2 costi 2 +2m 2 l 1 l 2 cos(8 1 -ti^ti^ (2) 

V = m 1 gl 1 cos8 1 + 2m 3 gl 1 cos8 1 + m 2 g(2l 1 costi 1 + l 2 cos8 2 ) (3) 

Thus the Lagrangianof the system is given 

L = -(m +m 1 +m 2 +m 3 )x 2 + (-m^ 2 + 2m 2 l 1 2 + 2m 3 l 1 2 ) Q 2 + - m 2 l 2 2 Q 2 + ("Mi +2m 2 l x + 
2m 3 l 1 )xti 1 costi 1 +m 2 l 2 x8 2 cos8 2 +2m 2 l 1 l 2 cos(8 1 -ti 2 )ti x e 2 - m 1 gl 1 cos8 1 - 2m 3 gl 1 cos8 1 - m 2 g(2l 1 cosd 1 + 
l 2 costi 2 ) (4) 

Diffcrentiatiim the La<rrandan bv and and yields Lasiraimc equation I ! } as: 

d dL dL _ 
dtd~8\~ dti[~ 
d dL dL _ 
d~td~8~ 2 ~ d~8~ 2 ~ 
Or explicitly: 

{^rn^ 2 +4m 2 l 1 2 + 4m 3 l 1 2 )d 1 + (m^ + 2m 2 l 1 + 2m 3 l 1 )x cos8 1 + 2m 2 l 1 l 2 8 2 cos(8 1 - 8 2 ) + 2m 2 l 1 l 2 8 2 2 sin(8 l - 
Qi) ~ ("Mi + 2m 2 /i + 2m 3 l 1 )g sm8 1 = (7) 

m 2 l 2 xcos8 2 +2m 2 l 1 l 2 8 1 cos(8 1 -8 2 )+^ m 2 l 2 2 8 2 - 2m 2 l 1 l 2 8 1 2 sin(8 1 - 8 2 ) - m 2 l 2 g sin8 2 =0 (8) 

Lagrange equation for the DICP system can be written in a more compact matrix form: 
D(9) 8+C(Q,8) 8 + G(9) = Hu (9) 

The stationary point of the system is (x, 8 1 ,8 2 ,x, 8 1 ,8 2 x) = (0, 0, 0, 0, 0, 0, 0), introduce small deviation around a stationary 
point and Taylor series expansion; the stable control process of the Double Inverted Pendulums are usually cos(9! - 9 2 )=1 



(5) 
(6) 
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)=Q[, sin9 2 S 2 . Linearization is made at balance position: \ 
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III. ANALYSIS OF DOUBLE INVERTED PENDULUM 

After obtaining the mathematical model ol the system features, we need to analyze the stability; controllability and 
Observability of systems in order to further understand the characteristics of the system [8]. 



A.Stability 

If the closed-loop poles are all located in the left half of "s" plane tin ; tern must b< tabl otherwise the system 
is unstable. In MATLAB, to strike a linear time-invariant system, the characteristic roots can be obtain by eig (a,b) function. 
According to the sufficient raid necessar} conditions for stability of the system, we can see the inverted pendulum system is 
unstable 



B. Controllability 

A system is said to be controllable if any initial state x(t ) or xO can be transfer to any final state x( t f ) in a 
time interval (t f - to), t > by some control u. 

The test of controllability due to Kalman if system is completely controllable if and only if the rank of the composite n 
Q c isn. 
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Q C =[BAB A nl B] (10) 

C. Observability 

A system is said to be observable if every state x can be exactly determined from the rr 
a finite interval of time < t > t f . 

The test of controllability due to Kalman if system is completely observable if and only if the rank of composite n 

O c =[C T A T C T (A T ) 2 C T (A T r'C T ] (11) 



"ulation. and tracking, 
indices [9]. 



IV. DESIGN OF LINEAR QUADRATIC REGULATOR 

This leads to the Linear Quadratic Regulator (LQRi system dealing with '-(ate regulation, output ix 
Broadly speaking, we are interested in the design of optimal linear systems with quadratic perfor 
We shall now consider the optimal regulator problem that, given the system equation. 
*(t)=Ax(t) +Bu(t) 
Y(t) =Cx(t) +Du(t) (12) 

Determine the matrix K of the optimal control h 
u(t) = -Kx(t) (13) 

So as to minimize the performance index 
J = £(x T (t)Qx(t) + u T (t)Ru(t))dt (14) 

Where Q is a positive-semi definite and R is a positive-definite matrix. The matrices Q and R determine the relative 
importance of the error. Here the elements of the matrix K are determined so as to minimize the performance index. 
Then u(t) = -Kx(t) = -R- 1 B T Px(t) is optimal for any initial x(0) state. 




Fig 2: Full state feedback representation of DIP 



Where P (t) is the solution of Riccati equation, K is the linear optimal feedback n 
Riccati equation. 

A T P + PA-PBR- 1 B T P + Q = (15) 

Where Q and R chose as Q = diag ([10 60 80 0]) and R = 1. 

Therefore, 

K= -fl-ifi 7 "? = [10 275.2453 -515.6502 16.2044 22.1046 -111.9285] 



e only need to solve the 



V. DESIGN OF PD CONTROLLER 

A Proportional-Derivative (PD) controller is a control loop feedback mechanism used in process control type 
industrial. A PD controller calculates an "error" value as the difference between a measured process variable and a desired 
set point. The controller attempts to minimize the error by adjusting the process input. PD controller calculation involves two 
parameters values proportional (P) and derivative (D). Proportional values is determines the reaction to the current error, 
derivative values determines the reaction rate at which the error has been changing [5~10]. 




Fig3: Block diagram of DIP PD controller 
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PD controller is combination of proportional plus derivative controller. It consists to a single input three output of 
a double inverted pendulum system. This is the case of stabilization of a double inverted pendulum at is x = r = 6 2 = 0, 
which is a physical unbalance position. So integral action will result in instability condition of control system and a simple 
PD controller is more adaptive. Then we have used a three (multi) PD controller. 

The transfer function of PD controller is (Kp + IQS) i.e. PD h PD 2 , and PD 3 transfer functions are respectively (Kj +K 2 S), 
(K 3 +K4S) and (K 5 +K 6 S). 




Fig4: The 



itrollcr 



We now use a pole placement technique of state feedback control system to determine the 6 PD control parameter. 
When assume that desired closed pole are -2.1+2.1425J, -5, -5, -5, -5. We can obtain the parameter of PD controller by using 
the MATLAB i.e. 

Ki=131.4941, K 2 =-23.4289, K 3 =-549.4167 

K4=-124.6183, K 5 =23.0785, K 6 =29.2327 



VI. SIMULATION AND RESULT 

The system under consideration and the proposed controllers are modeled and simulated in the 
MATLAB/Simulink environment. The step response performance of the two controllers is compared fig.5 shows the step 
response of the system under consideration in absence of a controller and is found to be unstable. 



Fig5: Step response of Double Inverted Pendulum without controller 
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Fig6: Step response of Pendulum angles and cart position by using LQR controller. 



Fig7. Step response of cart position by using PD controller 




Fig8: The response of first pendulum angle by PD controller 




Fig9: The response of second pendulum angle by PD controller 



The time response specification for the system under consideration equipped with the proposed controllers a 
given in Tables 2, 3 and 4. 
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TABLE 2 




Time response 


LQR 

controller 


PD 

controller 


Settling time 
(T s ) 


3.05 s 


3s 


Rise time (T r ) 


0.51s 


.15s 


Peak overshoot 


20% 


5.8% 


So State 


0.02 


0.0117 



F THE PERFORMANCE CHARACTERISTICS F< 





TABLE 3 




Time response 
specification 


LQR 

controller 


PD 

controller 


Settling time 
(T s ) 


4.68s 


2.876s 


Rise time (T r ) 


0.21s 


0.17s 


Peak overshoot 


1.6% 


7.5% 


Steady state 
error (ej 









CHARACTERISTICS Ft 



Summary of the 





TABLE 4 




Time response 


LQR 

controller 


PD 

controller 


Settling time (T s ) 


4.67s 


3.08s 


Rise time (T r ) 


0.907s 


0.62s 


Peak overshoot 


5% 


7% 


Steady state error 
(e ss ) 









clear that both are successfully designed but PD 



From both controller LQR and PD controller's re 
controller exhibits better response and performance. 

VII. CONCLUSION 

In this paper, LQR and PD controller are successfully designed for a Double Inverted Pendulum system. Based on 
he results, both controllers are capable of controlling the double inverted pendulum's angles and the cart position of the 
inearized system. However, the simulation result shows that PD controller has a better performance as compared to the LQR 
controller in controlling the Double Inverted Pendulum system. 

REFERENCES 

Yan xncli, JIANG hanhong. 'Control method study on a single inverted pendulum in a limulink environment". Measure and 

control technology, 2005, 24(07):37-39. 

ZHAN G Hongli "simulation of virtual inverted pendulum based on MATLAB". Machine Building and Automation, Dec 

2004,33(6): 103-105 

Chen Wei Ji, Fang Lei, and I i I i i i averted pc n rnd control in Aerospace 

and Astronautics, 2008. ISSCAA 2008. 2 nd international symposium on, vol. ,no.,pp.l -3,10-12 Dec 2008. 

Gugao Company, GT-400-SV inverted pendulum users ' guide. 2002. 

Yingjun Sang. Yuanyuan i an. Bin Liu. ' double inverted pendulum control based on three loop I'll) and improved 81' neural 

network second inlcrnalion intercut on digital niariulacmrin luioniation.201 1 ' 6 ' ' 

Bogdanov Alexander, 2004, "Optimal Control of a Double Inverted Pendulum on a Cart", Technical Report CSE-04-006.Y. 

A.L. Fradkov, P. Y. Guzenko, D. J. Hill, A. Y. 1 

Proc. of IFAC symposium on Control of Nonlinear Systems, Lake Tahoe. 1995:655-659. 

Jian Pan, Tun Wang "The study of two kinds control strategy based on inverted pendulum". Modern electronic 

technology.2008.1. 

Desineni subbaram Naidn. Optimal < ontrol system". Idaho state nni\ crsiiy. Pocatello. Idaho. USA, June 2002. 

Syed Hasan Saeed, "Automatic control system", S.K.KATARIA & SONS August 2002. 



